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ABSTRACT 

Solutions of type IIB supergravity which preserve half of the supersymme¬ 
tries have a dual description in terms of free fermions, as elucidated by the 
“bubbling AdS” construction of Lin, Lunin and Maldacena. In this paper we 
study the half-BPS geometry associated with a gas of free fermions in ther¬ 
modynamic equilibrium obeying the Fermi-Dirac distribution. We consider 
both regimes of low and high temperature. In the former case, we present a 
detailed computation of the ADM mass of the supergravity solution and hnd 
agreement with the thermal energy of the fermions. The solution has a naked 
null singularity and, by general arguments, is expected to develop a hnite area 
horizon once stringy corrections are included. By introducing a stretched hori¬ 
zon, we propose a way to match the entropy of the fermions with the entropy 
of the geometry in the low temperature regime. In the opposite limit of high 
temperature, the solution resembles a dilute gas of D3 branes. Also in this 
case the ADM mass of the geometry agrees with the thermal energy of the 
fermions. 
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1 Introduction 

According to the AdS/CFT correspondence, deformations of AdS geometries should map 
to states in the dual CFT living at the boundary of AdS [I] . Recently a concrete realization 
of this map has been found for the important sector of half-BPS operators of A/” = 4 Super 
Yang-Mills. These operators have conformal dimension A equal to the U{l)ji charge. They 
form a decoupled sector of A/” = 4 Super Yang-Mills which can be efficiently described 
by a gauged quantum mechanics matrix model with harmonic oscillator potential, see [ 2 ] 
and for earlier work [3]. The matrix model is well known to be completely integrable. 
The main reason behind integrability is that, in the eigenvalue basis, the eigenvalues 
behave as fermions in a harmonic potential. In the semiclassical limit the half-BPS states 
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can be depicted as droplets of fermions in a two-dimensional phase space. One expects 
then the following AdS/CFT dictionary. Small ripples above the Fermi sea correspond 
to graviton excitations of AdS^ x S^. Small holes below the Fermi energy correspond to 
giant gravitons, while small droplets of fermions ontside the Fermi sea map to dnal giant 
gravitons. All this is very reminiscent of both old and recent works on c = 1 string theory 
and its matrix model reinterpretation [1] [5] (for a recent review see |6]). 

Remarkably this whole picture has found an impressive confirmation through the ex¬ 
plicit construction of the full moduli space of half-BPS IIB supergravity solutions dis¬ 
covered by Lin, Lunin and Maldacena (LLM) [7] E. The phase space distribution of the 
matrix model eigenvalues is in one-to-one correspondence with IIB supergravity back¬ 
grounds which preserve half of the supersymmetry. Moreover the two-dimensional phase 
space of the fermions has an interesting physical embedding in the space-time geometry. 
At the quantum level the incompressibility of the droplets in phase space (due to Fermi- 
Dirac statistics) corresponds in the dual supergravity side to the requirement that the 
Ramond-Ramond five-form flux is quantized. The whole family of half-BPS geometries 
can be constructed in terms of an auxiliary function z which also determines the fermion 
distribution. The regularity of the supergravity background amounts to requiring a suit¬ 
able boundary condition on the auxiliary function. The AdS “bubbling geometries” are 
therefore in general smooth supergravity backgrounds. 

The fermions discussed so far are characterized by having a step-function distribution 
in the two-dimensional phase space. They can be seen therefore as fermions at zero “tem¬ 
perature”. It is then natural to investigate how turning on the temperature affects the 
supergravity solution. The fermion at non zero temperature are described by a Fermi- 
Dirac distribution. The corresponding AdS “bubbling” solution has been first obtained 
in [9] and further studied in [10] where it was given the name of hyperstar. This su¬ 
pergravity background can be thought of as resulting from a coarse graining process of 
smooth half-BPS geometries. The fermion distribution of the hyperstar fails to satisfy 
the boundary conditions necessary to obtain a smooth gravity solution. Quite generally 
when the smoothness condition is not satished naked singularities occur HB n. One 
expects that a' string corrections will modify the geometry in proximity of the singularity 
and that a horizon will be generated [T3| [13] . This class of singular supergravity solution 
can therefore be regarded as incipient black holes. An example of singular LLM solution 
is the superstar [15], which has been investigated from this point of view in Hg HB. 

^For related work see [8]. 
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The duality between fermion distributions and supergravity solutions at zero tempera¬ 
ture suggests that the thermodynamic properties of the fermion gas at hnite temperature 
should agree with the corresponding quantities in the supergravity side. In particular, one 
expects agreement between the thermal excitation energy of the fermions and the ADM 
mass of the supergravity solution. We will check that this is indeed the case in the two 
opposite regimes of low and high temperature. 

As we have already remarked, the hyperstar geometry is singular. The singularity is 
resolved quantum mechanically through the appearance of a hnite area horizon. One can 
then use the Bekenstein-Hawking formula to compute the associated entropy. By placing 
a stretched horizon in the hyperstar geometry we propose a way to match the supergravity 
entropy with the thermal entropy of the fermions in the low temperature regime, up to a 
numerical factor. 

Similarly we investigate the opposite regime of high temperature. In this limit the 
Fermi-Dirac distribution reduces to the classical Boltzmann distribution. We hnd that in 
this regime the metric is reminiscent of the so called dilute gas limit of LLM conhgurations 
associated to the Coulomb branch of Super Yang-Mills. 

The organization of the paper is the following. In the next section we present a 
concise description of the main results of LLM. In sec. 3 we review the thermodynamics 
of ID fermions in a harmonic potential. In sec. 4 we introduce the hyperstar geometry, 
discuss its properties in the low temperature regime, compute the ADM energy, angular 
momentum and entropy and we compare the results to the fermion prediction. In sec. 5 
we move on to consider the high temperature limit of the hyperstar distribution, we hnd 
the corresponding metric and determine its mass and angular momentum. We conclude 
discussing some open questions in sec. 6. 

2 Review of LLM 

In this paragraph we briehy review the LLM construction [7] of 1/2 BPS IIB supergravity 
backgrounds. These solutions correspond, in the dual CFT side, to states satisfying the 
BPS condition 

A = J (1) 

where A is the corresponding conformal dimension and J is a particular U{1) charge 
of the S'0(6) R-symmetry group. By selecting one generator of the S'0(6) R-symmetry 
group of A/” = 4 Super Yang-Mills we obtain a theory with S'0(4) x S'(9(4) x R bosonic 
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symmetry. In the dual supergravity description we look therefore for solutions with this 
isometry group. Assuming that axion and dilaton are constant and that only the self-dual 
hve-form held strength is turned on, the Ansatz for the background is 


= g^^dx^dx'" + ^dCl^ ( 2 ) 

F( 5 ) = F^ydx^ A dx'^ A dQ^ + F^^dx^ A dx'' A dfis (3) 


where the Greek indices /i, v run over 0,..., 3. The two three-spheres S'^ and in 
the metric make the S'0(4) x S'0(4) isometries manifest. The additional R isometry 
corresponds to the Hamiltonian A — J. 

For a background to be half-BPS there should exist a solution to the Killing spinor 
equation. Analyzing this equation, LLM were able to prove that the generic 1/2 BPS IIB 
supergravity background takes the form 


ds^ 

h-^ 

ydyV, 

z 

F 

F 

Bt 

dB 


—h ‘^{dt + Vidx^Y + h^{dy‘^ -|- dx^dx^) -|- ye^dVL\ -|- ye ^dVt\ 
2y coshG, 


€ijQjZ^ 


OjVj^ e^ijdyZ 


-tanh G 
2 

dBt A {dt -|- K) -|- BtdV + dB , 
dBt A {dt -|- K) -|- BtdV + dB 


1 3 ^ 

--y *^d 


z + 

y2 


Bt = 

~ 1 , 
dB = --y'^ *3 d 



(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


where i = 1, 2 and X 3 is the Hodge dual operator for the hat three-dimensional space 
parameterized by Xi,X 2 , y- Remarkably, the solution is completely specihed in terms of a 
single auxiliary function z{xi, X 2 ,y) which satishes the linear diherential equation 


didiZ + ydy 



0 . 


( 11 ) 


It is important to note that at y equals zero the product of the radii of the two privileged 
three-spheres is zero. Therefore, to avoid singular geometries, the auxiliary function 2 ; 
must satisfy a suitable boundary condition. This smoothness condition turns out to be 
= ±4 on the boundary plane y = t). In the limit z ^ 1/2 the sphere shrinks 
to zero while the other three-sphere remains hnite. The reverse statement applies when 
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2 ; ——1/2. It is conventional to assign black and white colors respectively to the 2 ; = —1/2 
and 2 ; = 1/2 points in the {xi,X 2 ) plane. If T> denotes a black region in this plane, the 
energy of the associated snpergravity solntion has the simple expression 

1 ^ ^ cPx^ ^ 


A = J = 


(Px 1 {xl + xl) 


'V 


2nh 


( 12 ) 


Jx) 2,Tih 2 h 

The plane has then a natnral interpretation as the phase space of one-dimensional 
fermions in a harmonic potential. This nicely matches the matrix model description in 
the dnal CFT side [2]. It emerges a beantifnl pictnre of the modnli space of half-BPS 
geometries of IIB snpergravity in terms of conhgurations of droplets of fermions on the 
(xi, X 2 ) plane. Note that the fnndamental equation ffTTj) has the symmetry z ^ —z which 
simply exchanges the S'^ and 5^ in the solution. In a held theory description of the 
fermions, this symmetry amounts to a particle-hole duality. 

The quantization condition on the total area A of the droplets is related to the hve- 
form hux N as follows 

with 


h = 27r/^. (14) 

The hux N coincides with the number of fermions. The simplest conhguration in phase 
space is a black circular droplet of radius Rq = \/2hN and the associated geometry is 
AdS^ X with N units of the hve-form hux. This background has A = J = 0 and 
corresponds to the fermion ground state. The boundary of the droplet can be thought 
of as the Fermi level of the fermions. The of the background is obtained by hbering 
the sphere on a two-dimensional surface S 2 in the {xi,X 2 ,y) space which encircles the 
droplet. One can easily obtain conhgurations with an arbitrary number of S'^’s by adding 
other droplets. If we deform the circular droplet to conhgurations with diherent shapes 
but same area, we obtain backgrounds with AdS^ x asymptotics. 

The fundamental equation flTTl) can be rewritten as a Laplace equation for the quantity 
^ = z/y^ in Si six-dimensional space with spherical symmetry in four of the coordinates. 
The coordinate y corresponds to the radial direction in the four-dimensional subspace. 
This observation reduces the task of hnding the full solution z{xi,X 2 ,y) of eq. (ITTll to a 
well known initial-value problem. Once the boundary condition z{xi,X2,0) on the y = 0 
plane is specihed, the solution is 


z{xi,X2,y) 



71 


z{x'^^ x'2i i^)dx'Px'2 

[(x — X ')2 -|- J/2]2 


(15) 
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We can similarly get 


Vi{xi,X2,y) = 


'-ij 

71 


z{x[, X 2 , 0 )( 2 :j — x'j)dx[dx 2 


[(x - x ')2 + j/ 2 ] 


212 


(16) 


Since we are going to consider only droplet configurations with radial symmetry, it 
will be convenient to rewrite the above formulas in polar coordinates {xi,X 2 ) —>■ {R,(j)). 
It is easy to see that in this case Vr = Vicos0 + V2sin0 = 0. Defining V = = 

R{—Vi sin0 + V 2 COS0) the differential equations relating z and V ([6]) read 


ydyV = -RdRZ, ^OrV = -dyZ. (17) 

R y 

Rewriting eq. (IT^ and eq. (ITbl) in polar coordinates yields 

z{R,y) = -1 z{R',0)-^zoiR,y;R')dR' (18) 


V{R,y) = j z{R!,Q)gv{R,y,R!)dR! 


(19) 


where 


Zo{R,y;R') 


gv{R,y]R!) 


- R'2 ^ y2 

2 [(R2 + R'2 + ^2)2 _ 4i?2^/2]l/2 

-2R^R'{R^ - R'^ + 1 / 2 ) 

[(R2 + R'2 + y2Y _ 4i?2^/2]3/2 ' 


( 20 ) 

( 21 ) 


We remark that Zq is the LLM function corresponding to a circular droplet. Indeed in this 
case z(R',0) = l/2Sign(i?' — Rq) and using eq. flTSl) one obtains z{R,y) = zo{R,y; Rq). 
As previously anticipated such a configuration gives rise to the AdS^ x solution. In 
fact performing the following change of coordinates [7] 


y = i?o sinh p sin 6* (22) 

R = Ro cosh p cos 0 (23) 

0 = 0 + f (24) 

one recovers the AdS^ x S'® metric in standard form 

ds‘^ = Ro[ — cosh^ pdt^ + dp^ + sinh^ pdDg + d9^ + cos^ Od^"^ + sin^ Od^l). (25) 


A variation of the method described so far can be similarly applied to obtain 1/2 BPS 
M-theory backgrounds with AdS 4 j x S*^’^ asymptotics [7]. In this case the geometry is 
in one-to-one correspondence with solutions of a three-dimensional Toda equation, which 
plays the same role as eq. mil. 
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3 ID fermions in the harmonic well 


In this section we review the basics of the thermodynamics of one-dimensional fermions 
in a harmonic potential. In what follows, we will consistently adopt nnits in which h = 
fcs = 1. We consider a gas of N non-interacting fermions with hamiltonian 

H{p,q) = ]^{p^+ q^) (26) 

in thermodynamic eqnilibrium at a given temperature T. For large N, we adopt the 
semi-classical approximation in which the energy is taken to be a continuous variable. 
The probability distribution as a function of the energy H{p, q) = e is given by the 
Fermi-Dirac distribution: 


where p is the Fermi energy. This is determined by the normalization condition 

f 

X efr-rt/r + 1 

which gives 


(27) 


(28) 


/^ = rin(e"''’’ - 1). 


(29) 


We will first consider the limit of very small temperature T, or more precisely N/T S> 1. 
In this limit, the Fermi level becomes 


p = N + 0{Te-^/^). 


The total energy of the Fermi gas is given by: 

E= ede 


Jo eh-M)/'^ -X 1 

which for small T can be evaluated by means of the Sommerfeld expansion 


(30) 


(31) 


/(e)* 


r-M 


TT 




1= I ■’ = / fU)de + —T^f(p) + - T^f'(p) + 0(T^). 

/o +1 Jo 6 360 ^ ^ 


This gives: 


E ~ 


7V2 


rjnZ 


(32) 


(33) 
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The first term is clearly the ground state energy of the N fermions, so we expect the 
dual gravity solution to have a mass (and angular momentum) difference of A = 
with respect to the AdS^ x background. It is worth noting that in eq. fl3^ we 
only neglect exponentially suppressed terms. In fact, since /(e) = e, there are no power 
series corrections to the energy beyond T^. This is a specific feature of the ID harmonic 
oscillator, and we will recover it in the energy § and angular momentum of the hyperstar 
in the low T limit. 

To evaluate the entropy of the fermion gas, it is convenient to first obtain the free- 
energy F of the system. This is computed from the partition function Z, which in the 
continuous limit we are considering reads 


Z 


exp 


f de ln(l + 

T Jo 


(34) 


One can verify that this expression for the partition function is correct by checking that 
the relation E = —^\nZ (where f3 = 1/T) is satisfied. Using the definition F = —TIn Z 
one obtains the free energy 


F 


Nfi-T de ln(l + 


N fi¬ 


de 


ed-AIT + 1 


= Nfi-E 


(35) 


where in the second line an integration by parts was made. The entropy is then given by 
the relation F = E — TS from which we get 




2E — Nfi 

T ■ 


For small T, using eq. and eq. one gets 

3 


(36) 


(37) 


where again only exponetially small terms are neglected. 

We now consider the opposite limit of very high temperature N/T -C 1. In this limit, 
the Fermi distribution clearly reduces to the Boltzmann density 


riFoie) nsie) = N(3 e 


^Modulo a subtlety involving terms to be discussed later. 


( 38 ) 








where (3 = 1/T. The total energy in this approximation is readily computed 

E = NT. (39) 

The entropy can be obtained from eq. (|36ll (which is valid for any temperature) using the 
large T approximation 

/i ~ T In N/T (40) 

and reads 

S N\nT + 2N - N\nN. (41) 

4 The hyperstar: Low temperature regime 

We now introduce the half-BPS geometry dual to the Fermi-Dirac gas described in the 
previous section [9]. This solution was named hyperstar in |10j . 

A given z{R, 0) corresponds to a fermion density n{R) in the phase space via the 
relation 

= (42) 

For example, the AdS^ x solution is associated to the step function density Uq = -d^R — 
Rq), which can be viewed as the zero temperature limit of the Fermi-Dirac distribution 
(EB- One can turn on the temperature on the fermion side by replacing uq with npoiR) 
and construct the corresponding supergravity background by using flTS]) . flT^ and fl42]) . It 
is important to remark that the temperature we are turning on is the temperature in the 
“auxiliary” description of the free fermion gas. It is not a temperature of the supergravity 
solution or of the dual gauge theory. Indeed, we remain in the supersymmetric half-BPS 
sector. It would be interesting to understand better what corresponds to this temperature 
on the gravity and gauge theory side. For the time being, we regard T just as a deformation 
parameter of the AdS^ x background. 

For low temperatures, the solution is a small perturbation of the circular droplet. In 
fact in this limit the fermion conhguration in the {xi,X 2 ) looks like a black disk with 
the boundary slightly “blurred”, as shown in £g. [H In the low temperature limit the 
expressions ([HD and ([inD can be obtained analytically as follows [9] 

1 ^oo o 

z'pDi.R^y) = 2^y nFDiR')-^zo{R,y;R')dR' 
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Figure 1: Droplet configuration in the low temperature limit. The greyscale ring around 
the Fermi level corresponds to a singular region of the spacetime. 


TT 


— Rq) +-^T 

0 


de^ 


ZQ{R,y]RQ = \Ple) 


e-Fo 
^ 2 


+ 0{T^) (43) 


and 


vId(R,v) = 


nFD{R')gv{R, y, R')dR' 


vr 


= Vo{R,y;Ro)--T^ 

b 


d (gv{R,y]Ro = V^e) 


de 




r 2 

P—1^0 
2 


+ C>(T^) (44) 


where we have used the Sommerfeld expansion fl32|) and where 


VoiR, y, Ro) — 2 


R^ + Rl + y^ 


[( i ?2 + i ?2 ^ j / 2)2 _ 4 ^ 2 ^ 2 ] 1/2 


(45) 


corresponds to the AdS^ x background. In these expressions Rq = \/2N is the radius 
of the droplet in the phase space at T = 0, and N is the number of fermions. It is easy 
to check that fl43l) and (1441) satisfy the differential equations flT7)l . 
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4.0.1 ADM form of the metric 


In order to compute the mass and the angular momentum associated to the hyperstar, it 
is convenient to perform the following change of coordinates 

„2 \ 1/2 


R = + 

y = Lr sin^ 

1 7 ^ 

4> = + 


COS 9 


(46) 


and we also rescale t ^ Lt to have conventional units. Here {t, r, Q^) parameterize the 
asymptotic AdS^ (in global coordinates), whereas (6*, 0, Hs) span the asymptotic S^. Of 
course, L is the radius of both the AdS^ and the S^. It is related to the radius Rq used 
by m via Rq = L^. In this system of coordinates the metric can be rewritten in ADM 
form as 

dr'^ / ~ .7 n 2 

- j 


ds^ = + gee 


+ dO"^^ + (d^ + M^dt) + gn^us (47) 


where Af is the lapse function and Af'^ is the shift vector. 
Introducing the expansion parameter 

27r^T^ tt^T^ 


7 = 


3L8 


6 A 2 


(48) 


and using the explicit expressions for and Vpp 


^FDidi,y) — 


R^-Rl + y^ 


2 [(i?2 + Rl + 1/2)2 - 4i?2i?2jl/2 

^ 2Ry{{Rl + y^ + R\Rl - y^) - 2R^) 

[(i?2 + RI+ 1/2)2 _ 4i?2i?2]5/2 


7 + 0 ( 7 ^ 


vfniR,y) = 7 


R^ + Rl + y^ 


[(i ?2 + Rl + 2 / 2)2 _ 4 ^ 2 ^ 2 ] 

? 2\2 


1/2 


1 + 


one obtains, upon implementing eq. 
here up to 0 ( 7 ^) terms § 

Af^ = il + 


2RiRmR^-Rir-yHR^ + Rl + 2y^)) 
l(ff + Bo + !/^)^ - iR^R^A'^ 

, the components of the metric, which we present 


L 2 


l-^L^F,{r,9) 


^We notice that our expression for differs from the one reported in [9]. 
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dee 


= 7 


2L{r^ + L^){r^ + cos‘^ 6) 
(r^ + sin^ 9)^ 

cos^e\l+^L^Fi{r,e) 


= 


gn^na - r 


^ 2 I r2 • 2fl -^2(r, 9 ) 

sin 9 


l-^L^F2ir,9) 




= L"sin"0 


l + -fL'^F2(r,ff) 


(50) 


where 

Fi{r,9) 

F2{r,9) 


(3 cos^ 9 — l)r'^ + 3L^ cos^ 9r‘^ + L^(2 cos^ 9 + sin^ 9) 

(r^ + sin^ 9)^ ’ 

(3 cos^ 9 — l)r^ + L^(3 cos"^ 9 — 2 sin^ 9)r'^ + L^(2 cos^ 9 — cos^ 9 — 1) 

(r^ + sin^ 9)^ 


A general property of LLM distributions with compact support is that the corresponding 
geometries are asymptotically AdS^ x S^. One can check that this remains true for the 
metric eq. (I50|) . This is consistent with the fact that the droplet of £g. [T]is effectively 
conhned in a hnite region of the phase space. From the expressions in fl5Tl) . we also notice 
that the Sommerfeld expansion seems no longer reliable in a region around the point 
r = 9 = 0. The appearance of eventual singularities will be discussed in the following 
section. 


4.0.2 Singularity 

The study of singularities for LLM geometries was undertaken in P and [12] @. There it 
was shown that all singularities appearing in the LLM supergravity solutions are naked 
and fall into two classes, namely timelike and null. While the former are considered highly 
pathological due to the presence of closed timelike curves, the latter are not. In fact, for 
half-BPS geometries with null singularity, the underlying fermion density function n{R) 

^The resolution of singularities in this context has been investigated in [10] . There, evidence was 
provided to support that spacetime singularities emerge from effectively integrating out the underlying 
quantum structure at the Planck scale. 
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always takes values in the region n{R) G (0,1). This is the case both for the hyperstar 
and the superstar solutions. 

To verify the presence of a singularity in the geometry, one should hnd a curvature 
invariant which diverges. The first non trivial invariant to consider is R\j^ = RmnR^^, 
{M,N = 0, ...,9), since the Ricci scalar R vanishes. Indeed, as a result of the Weyl 
invariance of the classical theory, the trace of the matter stress-tensor is identically zero. 
In order to check the consistency of the metric eq. (jhlljl we explicitly verified that this is 
the case. 

There are two ways to perform the computation of Rmj^- The direct approach involves 
the explicit form of the metric, while the indirect one makes use of the field equations of 
type IIB supergravity. In this case, the knowledge of the five form held strength ^(5) will 
suffice: 

Rmn = RNP1P2P3P4 ~ -^ 9 mnF‘^^ ■ (52) 

Suppose now we use the approximate solution for the metric, whose explicit form was 
given in the previous section, eq. fl50D . A lengthy calculation gives 


R 


2 _ 

MN — 


160 Pq^t, cos 9) 

+ siYi^ey 


+ 0{f). 


(53) 


In this 7 expansion, the hrst term corresponds to AdS^ x and is, of course, hnite. 
The linear term in 7 is, however, potentially divergent. Here PQ{r, cos 6) is a sixth order 
polynomial in both r and cos 9 and goes to zero as (r^ -|- sin^ 9) when r —>■ 0 and 
0 —0. The square of the Ricci tensor is therefore divergent when r = 0 and 9 = 0. 
It is interesting to note here that this is exactly the singular behavior one sees in Kerr 
black holes. Due to their angular momentum, the collapsing region is not a point but a 
zero-thickness ring. The Kretschmann invariant K RmnrsR^^^^ 1 for instance, for a 
Kerr black hole with mass M and angular momentum J = Ma, is 


^ 2 Qejr, cos 9) 

(r^ -f cos^ 9 )^ 


(54) 


where QQ{r, cos 9) also indicates a sixth order polynomial having the same behavior as 
PQ{r, cos 9) in the vicinity of r = 0 and 9 = 0. This is suggestive of the existence of an 
event horizon in the hyperstar geometry which may manifest itself through a'-corrections 
to the supergravity solution. 

This is not however the result one would have anticipated. From the form of the metric 
in LLM coordinates, it is quite natural to expect a singularity at y = 0. Using flT 6 |l . we 
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can see that this corresponds to r = 0 or 0 = 0 in asymptotic AdS^ x coordinates. On 
the other hand, the singnlar region appearing in fl53l) is mapped to {R = L,y = 0), which 
is jnst the Fermi snrface of the fermions. We expect the singularity to be at least smeared 
over an extended region around the Fermi energy, since there the fermion density is less 
than one, see £g. [H 

What is therefore the true singular region of the hyperstar? We can try to address 
this question in a quite general fashion valid for all LLM geometries. We simply need 
to know the behavior of the functions z{R, y) and V{R, y) in proximity of ?/ = 0. We 
can distinguish two different cases depending on whether zo(-R)O) = \miy^Qz{R,y) is 
independent of the radial coordinate R or not. In what follows we will focus on the latter, 
since this is the case of the hyperstar. 

We would like to hnd z{R, y) and V(R, y) in terms of an expansion in y or functions 
of ?/, such that the differential equations ffT7|l will be order by order satished. It turns out 
that the appropriate Ansatz is the following 

z{R,y) = zo{R,Q) + fi{R) y'^ \ny + ... 

V{R,y) = Vo{R,0) + gi{R) \ny + ... (55) 

The functions fi{R) and gi{R) are determined to each order from the same differential 
equation (fT7|) . For the case that concerns us here we have fi{R) = —^dR^RdnZo^R, 0 )), 
Vo(i?, 0) = —^RdjiZo and gi{R) = —ROrZq. It is now easy to hnd the complete solution 
for the metric and the hve-form held strength in this region and subsequently calculate 
^MNi using either of the methods indicated above. We hnd 

^MN = ——^2-^ ^2 (^o(-R, 0 ), fi{R)) Iny + ... (56) 

y 

where hi(zo(-R, 0)) and h 2 {zo{R, 0), fi{R)) are non-zero functions of the variables indi¬ 
cated. 

Indeed we see that the leading term is divergent ai y = 0, as expected. We must 
therefore conclude that this is the singular region of the hyperstar, and that we cannot 
rely on the Sommerfeld expansion fl3^ for calculations in the small y region. This will be 
important later for computing the entropy through the Bekenstein-Hawking formula. 

4.0.3 Flux 

To check the consistency of the hyperstar solution, we can verify that the hux of the 
hve form ^( 5 ) remains equal to N, independently from the temperature T of the fermion 
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gas. From general considerations this has to be expected, since the temperatnre can be 
viewed as a tnnable continuous parameter and as such it cannot modify the flux which 
is a topological constraint. At zero temperature, i.e. for the AdS^ x solution, using 
the explicit expressions for the held strength in the LLM solution and the change of 
coordinates eq. (l46jl one obtains 

3 

= — dt A dr A dQ^ + 2N sin^ 9 cos 9 d9 A d^ A dQs. (57) 

The hux is computed by integrating over the S'®, and including the appropriate 

normalization is equal to N. To check that temperature perturbations do not alter the 
hux, one has to verify that corrections to vanish when integrated over the hve 

sphere. Up to second order in the temperature, these take the form 


-- 

{5)9<f>Q3 

^( 2 ) 

(5)e^n3 


■ 3n n L^P6{r,COs9) 

7 sm 9 cos 9 ——---— 

(r^ + sin^ 9 )"^ 

2 • 3/1 /I L^^pio{r,cos9) 

7 sm 9 cos 9 ——---— 

(r^ + sin^ 6 *)® 


(58) 


where pe and pio are polynomials of degree 6 and 10 respectively. Although the explicit 
formulas look rather involved, the integration over 9 can be carried out exactly at arbitrary 
r and indeed yields 








0 . 


(59) 


4.1 Mass 


In this section we present a systematic derivation of the ADM mass of the hyperstar 
solution. The natural expectation is that this mass should coincide with the thermal 
energy of the auxiliary fermion gas system. 

The Einstein-Hilbert action in a d-dimensional spacetime is 


Sgrav = {R - 2A) - (f d'^ (60) 

IbirGd Jm oirGd Jom 

where we included the Gibbons-Hawking boundary term and is the metric on the 
{d — l)-dimensional timelike boundary. Following [19], the quasi-local stress-tensor can 
be computed by the variation of the gravitational action with respect to the boundary 
metric 




2 6S. 


grav 


-7 <^ 7 , 




(61) 
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Using (1601) this is 


^ - ©7^'^). (62) 
SvrGrf 

In the previons expression we have introdnced the extrinsic curvature of the {d — 1)- 
dimensional timelike boundary embedded in A4 

e"" = (63) 

and we denoted the corresponding trace by ©. The covariant derivative is taken with 
respect to the metric g^y of the full spacetime and if is the unit normal to the boundary. 
The stress-tensor fl6T|) generically diverges as we approach the boundary dM. when the 
spacetime is asymptotically AdS. In the context of the AdS/CFT correspondence we can 
view the gravitational quasi-local stress-tensor as the expectation value of the stress-tensor 
in the associated conformal held theory. The divergences get then a natural interpretation 
as standard ultraviolet divergences in quantum held theory [20]. We can regularize the 
theory by adding suitable counterterms to the original stress-tensor 

i • (64) 

87 rGd V f-'yhiiiyj 

The counterterms are consistently constructed using only the boundary metric j^y and 
its covariant derivatives and are (almost) uniquely determined by requiring a cancellation 
of the divergences and general covariance (for a review, see [21]). The boundary metric 
'jfj^y can be written in the ADM form 

'jfj^ydx^dx^ = —M^dt^ + aabidx°‘ + J\f^dt){dx’^ + M^dt) (65) 


where S is a surface of constant t inside dJfl. Conserved charges are obtained by inte¬ 
grating over a spacelike hypersurface at inhnity. A hnite expression for the mass is 
obtained substituting the regularized stress-energy tensor in the following formula 


M= d^-’^xs/aM^.u^T^yv!', 


( 66 ) 


where is the timelike unit normal to S. For instance, the application of this method 
to the hve-dimensional AdS-Schwarzschild black hole 


ds^ = — 


_ I 1 _ 

L2 ^ V r 


dr + T 




ih. -L 1 _ f 10 V 

L2 \r ) 


j+ r‘^{d9‘^ + sin^ 9dcf)"^ + cos^ 9d'ljj'^) (67) 
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yields [ 20 ] 


Svr/^ 


Ovrrg 


M = ( 68 ) 

32G5 SGs 

The first term, which is present also when the black hole disappears, corresponds to the 
Casimir energy of the vacuum in the dual CFT. 

It would be nice to have a similar counterterm method directly in a ten-dimensional 
setting. Unfortunately, extending the program of holographic renormalization to ten¬ 
dimensional metrics with AdS^ x asymptotics seems problematic [22] . We are therefore 
forced to use alternative approaches. In the first one, we will determine the relevant 
components of the stress-tensor relative to some reference geometry following [23]. The 
second approach is the so called background subtraction method [23]. In both cases one 
has to carefully match the asymptotic geometry of the supergravity solution with that of 
a reference background. Neither of the methods can reproduce the Casimir energy of the 
associated CFT. However this will not be a problem in our case since we are interested 
in computing the energy difference between the half-BPS supergravity solution and the 
AdS^ X ground state. 

We now proceed to compute the mass of the hyperstar (1501) as a series expansion in 
the small parameter 7 = This mass should agree with the energy of the free fermion 
gas, eq. fl3^ . We will first consider the leading order in 7 and comment on 7 ^ orders in 
a later section. 


4.1.1 First approach 

Following [23], we obtain the stress-tensor associated with the metric (1301) relative to the 
AdS^ X background metric We need to require that the difference between the 
two metrics falls off suitably fast for large radius. Explicitly we want that 

Qrr-glr = o(l/r®) 

9ra-gla = o{l/r^) (69) 

where o(l/r"') means that these differences go to zero more rapidly than 1/r" and the 
index a runs over all the coordinates except r. To satisfy such requirement we implement 
an appropriate change of coordinates (r, 6 ^) —>• (r, 6 ^), which we presently discuss. The 
effect of using these new coordinates is to make the leading asymptotic perturbations of 
the metric all in components parallel to the boundary directions. Then the line element 
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becomes 


ds'^ = + '^^dx'^dx^ + ... ( 70 ) 

from which one can read off the stress-tensor up to a multiplicative constant depending 
only on the space-time dimensions. 

The first step is therefore to find a coordinate system such that the metric satisfies 
We consider the Ansatz 


r 

cos 9 


~ fl f2 
H—z—1“ 

n I fs I 

cos9 + — + — 


(71) 


where the fi = fi{9) (i = 1,... ,4) are functions to be determined in order to adjust the 
asymptotics of the metric. 

In terms of the new variables r and 6, the Qff component of the metric has an expansion 
for large r which differs from the background reference metric g% = {l + r/L‘^)~^ by terms 
containing the /j. The 1/r'^ term can be eliminated by tuning /i = |(3cos6*^ — 1)T^7 
and similarly the 1/r® with an appropriate choice of f 2 - The first constraint in (l6^ is 
then satisfied. Analogously, /s and are fixed by requiring the vanishing of the l/r^ and 
1/r^ terms in g^g, which appears after changing variables according to flTTl) . Once the fi 
are fixed, one can verify that the other components of the metric coincide with g^^ up to 
orders O (7“^). The only disagreement is found in gu and guifii, which contain a term at 
order O (r°) 


7(3 cos^ 6^ — 1) 


( 72 ) 


which, nonetheless, vanishes upon integration over the (including the appropriate 
measure). We notice that the same factor already appeared at leading order in the 
asymptotic expansion of the metric perturbation, see eq. 0511) . 

From 07Up and the explicit expression for 

= + ( 73 ) 

one can read off the time-time component of the stress-tensor 

Ttt = + l + 


18 


(74) 


^ ^4(3 cos^ 9 — l)r^ + L^(ll — 39 cos^ 9 + 60 cos^ 9)^+0 ■ 


This expression has to be integrated at the spacelike boundary in order to give the mass 


M 


4 

IOttGio 

4 

IbvrGio 


j + Ttt 
L^2tt){2ttY 


t/2 


d9 cos 9 sin^ 9 Tf 


tt 


(75) 


where Gio = and fi = r ^ y^fi '00 9oe gfm is the integration measure. The hnal 
result for the mass is 


M = 




(76) 


which agrees with the thermal excitation energy of the N fermions above the ground 
state, eq. ([33]). The extra L in the denominator comes from the rescaling of the time 
variable already discussed. It is important to remark that in obtaining these expressions 
we have consistently worked at order 7 . We will comment on the signihcance of higher 
order terms in a later section. 


4.1.2 The superstar 


As a further check of the validity of the procedure just discussed, we also apply it to the 
so-called superstar, a family of asymptotically AdS^ x backgrounds discovered in [15] 
and further studied from the LLM perspective in [16] [TT] [ 1 ^ [10]. The extremal 1/2 
BPS superstar metric is governed by two parameters, the flux N of the 5-form through 
the and one of the three angular momenta on the S^, J 3 , which coincides with the 
energy A because of the BPS condition. Explicitly the metric can be written as HD 


ds^ 




2 \ T^TJ T 

^ -I——— sin^ 9d(j)'^ + 9dtd(f)+ 

(j (j 


+G ( -I- r‘^d+t^ + L‘^Gd9^ -|- ^ cos^ ddflg 


(77) 


with 


f = l + H 


L 2 ’ 


G = \/ sin^ 9 + H cos^ 9 , H = 


l + + (78) 
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Also in this example we want to satisfy the fall off conditions fl69ll . By choosing an 
appropriate coordinate system as in flTTl) it is easy to see that 


ftt 


^(2-3cos2 0)V+ 

+ 6^2 ~ — (4 + 15 cos^ 0 — 24 cos"^ 0)Q) + (9 j . (79) 


The expression for the mass is then 


M 


4 

IOttGio 


f^Ttt 


4 

IhvrGio 




A 

L 


(80) 


which, up to the L coming from the rescaling of the time, is exactly the energy of the 
geometry. Note that we have again neglected contributions quadratic in Q in the stress- 
tensor fl79l) . 


4.1.3 Second approach: Background subtraction 

We now discuss the second approach |23] for computing the mass of the hyperstar. In 
the background subtraction prescription the ADM mass is obtained by integrating the 
quasi-local energy Afi^K — Kq) over the (d — 2)-dimensional spacelike hypersurface S at 
radial inhnity 

To obtain M one needs 77^^, the extrinsic curvature of S embedded in a constant time 
hypersurface 

A'-- = (82) 

Now the covariant derivative is calculated with respect to the metric of the constant 
time hypersurface, and = grr ^ <5^ • In flSTl) K and Kq are the traces of the extrinsic 
curvature of the spacetime and of the reference background respectively, and /i is the 
measure on S. 

In this case we also need to carefully tune the components of the boundary metric 
with those of the AdS^ x background by performing an asymptotic coordinate trans¬ 
formation as in the Ansatz dZU). Let us consider the extremal superstar solution in its 
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five-dimensional reduction to understand which fall-off requirements we need to impose. 
The mass of this solution was hrst obtained in [25j. The line element reads 

ds‘^ = {f~^dr^ + r^dQl) (83) 


where H and / are defined as in (|78|l . The parameter Q appearing in (178il is the five¬ 
dimensional electric charge and corresponds to the angular momentum J in the ten¬ 
dimensional uplifting of the superstar solution fITTD . see also [26]. We perform the following 


change of variable on the solution 

r~2 = 

(84) 

which asymptotically amounts to 

r = f-^ 

6f 

(85) 


A posteriori one can verify that additional higher order terms in (l8^ do not modify 
the final answer for the mass. After this transformation, the difference between the 
components of the boundary metric and the global AdS^ background becomes of 

order (9(r“^). An explicit calculation of the extrinsic curvature yields 


K = 


3 3L 
L~^ 






( 86 ) 


To obtain a hnite mass, we need to subtract the extrinsic curvature of AdS^ 


Ko = -- 1 + 


2\ 1/2 


(87) 


Using 


M = 


SvrGn 


I^M{K-Ko) 


( 88 ) 


r/4i 


with Gs = 7r/4 □, we obtain the well known result 


M = Q. 


(89) 


One can easily verify that if we had not implemented the transformation fl85l) we would 
have gotten 


K = 


3 _ 3L 
L~^ 


3L3 Q2 3^gx 

~8~ si ~Y~) 



( 90 ) 


®In this example we use the units of [25] . 
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and correspondingly the incorrect result 


M = IQ- ( 91 ) 

As in the ten-dimensional example fITTD . we have again neglected a term proportional to 

Q2. 

We now proceed similarly with the hyperstar solution using /i, /2 to £x the asymptotic 
behavior of and analogously /a, /4 to fix B Having four parameters at our 

disposal we require that and are of order O ■ With this choice we 

obtain Sgge = 0{f~^), while for the other component of the boundary metric g^^^ we have 
6g^^ = 7L'^(—1 -|- 3cos(6')) which integrates to zero on the S^. After having implemented 
this coordinate transformation, we can compute the extrinsic curvature to linear order in 
7 obtaining 

^ ^ -|-^(3-7(3cos20-1)7) + 

+ ^(3 + 4(28 - 159cos2 0 + 174cos^0)7) + C> . (92) 


Subtracting the extrinsic curvature contribution of the background 


Ko = -- 1 + 


2\ 1/2 


3 3L 3L^ ( 1 

+ —+ C> 


L \ ) L 2r‘^ 8r‘^ 

and using the ADM mass formula, eq. (l8T|l . we obtain 

M = 

4 ' 6L 


(93) 


(94) 


which is again the expected result. 


4.1.4 Contributions to the mass of order 7 ^ 

It remains to discuss the relevance of the quadratic terms in 7 that we have so far con¬ 
sistently neglected. According to the discussion following eq. (I33|) . we would not expect 

®We could have also chosen to use the parameters ft to fix the other components 5^0 ,ggg of the 
boundary metric. This ambiguity alters only the quadratic contribution to the mass which, as will be 
discussed, is not physical. 
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contributions to the mass at orders higher than j T‘^. We now check whether this is 
the case. Using the expressions at order 7^ ~ for Zpjj and Vp^ 


zlD\R^y) 


U, 


T (2) 


FD 


{R,y) 


^2__ 

^ 5 [(/?2 + Rl + J/2)2 _ AR^RlfR 
{{Rl + + R\Rl - lly^){Rl + y^f+ 

+ mly^^ - 2y^) + R^{llRl + lAy"^) - 

2 84i?8 i?2 

^ 5 [(/?2 + RI + J/2)2 _ 4/?2i?2]9/2 

((i?2 - 4j/2)(i?2 + y2f _ ^ 9^4^2 _ 9^2^4 _ ^4^6)^ 

+R\<oRl + 21i?2|/2 + 6^4^ _ ^6^4^2 ^ ^^^2^ ^ ^95) 


it is straightforward to write down the corresponding asymptotic expression for large r 
of the hyperstar metric, which is not particularly illuminating and, therefore, we do not 
present it. 

It is not difficult to see that, differently from what expected, there seems to be a non 
vanishing contribution to the mass proportional to 7^ ll. The exact coefficient of this term 
depends on the procedure used to compute it. In the hrst approach discussed above there 
is a quadratic contribution to the stress-tensor 

2 

f ® = -^ (19 - 159 cos^ e + 216 cos^ 9) (96) 

and to the mass 




2 


The method of background subtraction gives 


TT 


72L9 




(97) 


and 




540 cos^ 9 + 747 cos"^ 9) 


(98) 




125L^ 2 ^ 1257r^^4 

128 ^ 288L9 


(99) 


^On the other hand, the angular momentum, which can be obtained from ~ ^ does 

not receive corrections beyond 0 ( 7 ). 


23 









The presence of this term and its scheme dependence are, however, not completely 
surprising and have already been discussed in the literature. In computing the superstar 
mass, both in hve and ten dimensions, we already encountered a similar issue, see eqs. 
fl7^ and fl86l) . Indeed, retaining the terms in the computation of the mass, one would 
obtain a non linear BPS condition M Q — ^ This relation clearly conflicts with 
the expectation M > |Q|. One can nevertheless recover the usual linear BPS condition 
by including appropriate finite counterterms related to scalar helds [2H]. This discussion 
can be generalized to the three-charged AdS^ black hole. It has been observed in [29] that 
terms quadratic in the charges are related to a trace anomaly of the stress-tensor. This 
anomaly stems from a renormalization scheme which violates the asymptotic isometry 
group of AdS^ and can be removed by adding to the action the hnite counterterm proposed 
in [28|. 

In the light of these examples, we therefore consider fl97|) and fl9^ as spurious: They 
should be eliminated by a convenient choice of counterterms, although we do not know 
how to carry out this procedure directly in ten dimensions. 

Orders beyond 7^ do not contribute to the mass of the solution, because they fall off 
too fast at radial inhnity. 


4.2 Angular momentum 


As a check of the BPS condition for the hyperstar solution, we now calculate the associated 
angular momentum. This computation is most easily done in a five-dimensional setting. 
The ten-dimensional angular momentum J coincides with the electric charge Q of the 
U{1) gauge held A coming from dimensional reduction on the S^. The gauge held can be 
read oh from the term g;j^^{d(j) + M^dtY in the ADM metric and therefore coincides with 
the shift vector: 


2L 

A = dt = —7 dt + O 

The associated charge (angular momentum) is then 

J = - / *5 dA = 7— = — 

IbTrGs J S3 4 6 


( 100 ) 


( 101 ) 


where *5 is the hve-dimensional Hodge star operator. In our normalization the hve- 
dimensional Newton constant is G5 = Giq/VoI{S^) = The LA factor in fllOip is 

^Looking at the five-dimensional gauged supergravity action one would expect a contribution to the 
charge of the type A A F. This term is nonetheless subleading and vanishes at radial infinity. 
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necessary for obtaining conventional units. Comparing J with the mass formula eq. fl76l) . 
we obtain the BPS relation M = J jL. 

4.3 Entropy 

In the previous sections we have found agreement between the ADM mass and the ther¬ 
mal energy of the fermions. Since the Fermi gas has non-vanishing entropy at non-zero 
temperature, we expect the same to occur for the supergravity solution. We would like to 
understand how this entropy arises geometrically in the case of the hyperstar. Although 
the solution we are considering seems to have a naked singularity, it is expected that 
a' corrections to the equations of motion might generate a hnite-area stretched horizon. 
With these corrections we can think of the hyperstar as a legitimate black hole. 

In the presence of an event horizon, the entropy of a gravitational solution in d di¬ 
mensions is given by the celebrated Bekenstein-Hawking formula 

5 =^ ( 102 ) 

where Ad is the area of the horizon. In our case the entropy is still given by fll02p but 
now A = Ash is the area of the stretched horizon. 

Since we do not know the explicit form of the a' corrections, the location of the 
stretched horizon is inherently ambiguous. Therefore we expect to reproduce the fermion 
entropy up to a numerical coefficient. 

As we already discussed, the y = 0 plane is a null singular region. It is reasonable to 
assume that the a' corrections will generate a horizon at ygh — 0 -|- 0{a'). 

We therefore need to compute the area of the y = ysh plane, with ysh — a' = gs Ip ~ 
gs in units where h = 1. This area turns out to be hnite. The metric in LLM 
coordinates for hxed t and y reads 

ds^\t^y = fixed = dcf? + h‘^{dR‘^ + R^dcfP') + ye^dVL\ + ye~^d^l ( 103 ) 

so that the integration measure is 

p = Ah^R^ -V^y^zz h^Ry^ “ (i “ ^ 

where we have assumed the expansion 0551) . so that the term V'^ ~ In^?/ can be neglected 
for small y against ~ y~^ and ~ zq. By restricting the measure 01041) to y = ysh, the 
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Bekenstein-Hawking formula yields 


^ = 


^sh{y = Vsh) 27r(27r2) 


AG 


10 


4-2714 


2\2 roo /t 

ylh RdR 


1/2 


~ c / RdR\/n{l — n) 

Jo 


( 105 ) 


where c is a numerical constant. For the hyperstar n = upn so that 

ynpoil - npo) = 

with e = R^/2 and /i = T — l). Using eq. fllOSp we obtain 

r°° 

Sc^c de- - 

Jo 1 + 

/TT _iR\ / TT 1 

= 2cT — — arctane 2 =2cT — — arctan , ^ . . 

V2 J V2 Ve^/T_i; 

In the low temperature approximation this yield 

S oc T (1 + 0(6-^'’')) . 


(106) 


(107) 


(108) 


Therefore the entropy is proportional to T, as expected from eq. (ETj), up to corrections 
which are exponentially suppressed for 3> T. 

In the high temperature limit, however, eq. fll07p does not seem to reproduce the 
logarithmic behavior of the Boltzmann entropy. In this limit, which will be studied in the 
next section, the assumptions and the approximations which led to eq. fllOSp might not 
be valid since T is not a small parameter. 


5 High temperature regime 

We now move to consider the high temperature regime. In this limit the Fermi-Dirac 
distribution reduces to the classical Boltzmann distribution. Correspondingly, the droplet 
spreads over a larger part of the y = 0 plane and the singular greyscale region is not 
conhned inside a thin ring anymore, as shown in £g. [2] . 

The auxiliary function z can be computed in this regime as 

17“ d 

^liR^y) = 2^ J ^BiR')-^zoiR,y,R')dR' (109) 
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Figure 2; Droplet configuration in the high temperature limit. 


where 

nsiR) = Nj3e~^^ 


( 110 ) 


and (3 = 1/T. Making the change of variable = e and using the explicit expression 

for zo{R,y, R!) we can write the integral as 


zl{R,y) = --2y^NPh{P,R,y) 


where we have dehned 


R{j3,R,y) = dee 


(2e + y^ + R^) 


[(2e + |/2 + i?2)2_3^2g]3/2' 


( 111 ) 


( 112 ) 


The high temperature limit corresponds to the small (3 region. Therefore we want to 
hnd an approximate expression for Iz{j3,R,y) near (3 = 0. It is easy to verify that 
Iz{0,R,y) = and also that 


diz _ r ee-^\2e + R^ + y^) 

Jo [(2e + 1/2 + R 2)2 _ 8i?2g]3/2 


(113) 


diverges as In (3 in proximity of (3 


0, because the integrand goes like e je for large e. 
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This suggests a low (3 expansion of the form 


= R{0,R,y) + Aplnp + ... 
zl{R,y) = ^-2y‘^Nf3 (^-^ + Af3\ni3+ ...) . 


(114) 


Since it is not possible to compute explicitly ^ for /3 —0 because of the divergence, to 
hnd its small (3 behavior we find it useful to hrst regulate the integral by considering the 
quantity 


dh , 1 


e e 


-/3e 


e2 + (i?2 + y2y 


The new piece 


In = -- 


de 


e e 


-3^ 


e2 + (i?2 + y2y 


(115) 


(116) 


has the same divergence structure of ^ and its value is known for finite (3 in terms of the 
Sine and Cosine Integral functions Si{x), Ci{x). The corresponding small (3 expansion 
can be given explicitely as 


+ \hi3(k!‘ + v3] - \l3(R" + a") + O(tsH^ts) 


(117) 


where 7 is the Euler-Mascheroni constant. The combination (I115p is by construction 
convergent for any /3, and has a well defined /3 —> 0 limit which can be easily computed 
analytically 


diz 

d(3 


-h 




/ 3=0 


+ bny^ - iln(fl2 +j/=). 


(118) 


Using eq. flll7p we obtain the high temperature expansion of ^ 


+ In 2/2 +0(/31n/3) 


^ = \\n(3 + \U-\n2 + l-^ 

5/3 4 4 V y 

and integrating in (3 we can finally get 

W.R,!/) = ^ + jTJlii/S + ^ It - ln2 - + ln!,q + C>(y 111/3). 


(119) 


( 120 ) 


The corresponding high temperature limit of z]^, keeping only the first two orders, reads 
then as follows 

1 


zM^,y) = --N/3 


Ny^ 


(3^ ln/3 + 0(/32). 


( 121 ) 











To obtain the metric we need to find also the fnnction V^{R,y). Starting from eq. 
(IT^ and inserting the Boltzmann distribution we arrive at 

2R'^{2e + + R^) 


V^{R,y)= deN(5e 


-Pe_ 


[(2e + 1/2 + i?2)2 _ 8i?2g]3/2 


= 2N(3R^Iv{(3,R,y). (122) 


One can verify that Iv{0,R,y) vanishes and that diverges logarithmically in the 
(3^0 limit. We can proceed similarly as before by regulating with an appropriate 
“reference” integral, to finally obtain 


Iv{/3,R,y) = ~p\n/3-^{^-\n2 + l + W) + OipHn/3). 
In proximity of /3 = 0 the leading contribution to is therefore 

VS(R,y) = -l]VflVln/3+ 0(/3"). 


(123) 


(124) 


The expressions for and Vj consistently satisfy eq. flTTj) . Note that does not depend 
on R and that similarly does not depend on y. This fact is nevertheless an artefact 
of the approximation we made. To study its limits of validity, we can look at (I12ip and 
(I124j) and require that the corrections are small: From this we can infer the conditions 

1 1 




R^ < 


/31n/3’ “ "" Nf3nnf3' 

We can now find the metric at first order in the low (3 expansion, i.e. = 1/2 — Nj3 
and Vj = 0. The metric in the LLM coordinates is quickly computed and reads 


ds^ = 


y 


i^—dt^ + dflg) + 


y 


VW 

Rescaling the coordinates as 

i = V = 

brings the metric into the form 

1 1 


+ y^d^l + + R^d^^ 


R = {Npy/^R 


(126) 


(127) 


ds^ = {NfdY^^y ^—dt^ + 


jdflg I + 


+ fdnl + dR^ + R^d^‘^) . (128) 


' {N(3y/^y 

This form of the metric closely resembles the dilute gas approximation limit studied in 
[7]. There, one considers a configuration of droplets with area Ai in the {xi,X 2 ) plane, 
and send the distance between the droplets to infinity by the rescaling 


X 


Xx, 


X 


Xx, 


y ^ >^y, A 


cxo 


(129) 
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while keeping the droplets areas Ai hxed. The corresponding metric reads 


ds^ = H-^/^ [-dP + x^dnl] + 

where the harmonic function H is 


+ y^dVti + dx^dx^ 




A 


^ I l(i - x'X + y^P' 


(130) 


(131) 


Thus the metric eq. flldOp can be viewed as a multi-center solution for a stack of separated 
DS-branes, and corresponds to the SO {4) invariant sector of the Coulomb branch of the 
gauge theory. 

Upon the identihcation A = one can see that the dilute gas limit A —> cx) 

is similar to the high temperature regime f3 —>■ 0 of the thermal solution eq. (11281) . This 
is perhaps not surprising since in the high temperature limit the fermion density goes to 
zero. We also notice that a continuum version of eq. fll3ip with = d'^x'/yJNfd gives 


1 j' d'^x' 1 1 

Try ^/N(3 [(£ - x'Y + y/NPy^ 


(132) 


which is what we would expect in order to match eq. (11281) with eq. fll30p . 

Taking into account the next to leading order corrections for z]^ and in eq. fll2ip 
and (11241) . we obtain the metric 


ds^ = \-dt^ + (1 - Np)dnl \ + [dy^ + (1 + N/3)y^dnl + dR^ + R^dcj)'^] 

+^R^yp'^^‘^\nj3dtd(t) (133) 


where 


H 


N(5 

y2 


iV2/32 

y2 


+ l/3Wln/3. 


(134) 


At this order we have a non vanishing and this determines the presence of the mixed 
term g^t in the metric. 


5.1 Energy and angular momentum 

We remark that the region of validity of the approximations made so far does not allow 
us to use the metrics (I126p and (I133p in the asymptotic region R^ + y"^ 1, because of 

the conditions eq. (I125p . Therefore, to compute the energy of the hyperstar in the high 
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temperature regime, we need to find the form of the metric in the complementary region 
of validity. The new metric will be trustable in the asymptotic region and will allow a 
calculation of the energy with the methods already discussed. To this end, it is convenient 
to hrst introduce polar coordinates in the {xi,X 2 iy) space 


R = u cos'd 

y = Msin'd. (135) 


Then one can evaluate eq. (Illip and eq. fll22p in an expansion for u S> 1 while keeping 
T hxed but large (such that we are in the Boltzmann regime). The integrals involved in 
the expansion can be readily computed analytically and one ends up with the result 

z^iu^d) = - —2iVsin^— 8 iVT sin^ ^(3 cos^ d — 1)^ 

2 

—48iVT^sin^'d( 10cost'd — 8cost'd + 1)^ + O (^ ) (136) 

\u^ J 


Vj(n,d) 


2N cos^ + 8 NT cos^ ^(3 cos^ — 2) \ 

+48 iVT^ cos^ ^(10 cos"^ -d — 12 cos^ d + 3)^ + C> 



(137) 


where in the expansion we have kept only terms which contribute to the mass and angular 
momentum. One can now go to the AdS^ x coordinates via the change of variables 
given in eq. fl46p and use z'^ and to obtain the asymptotic form of the metric. The 
explicit expressions are somewhat lengthy and we will not report them here in detail. 
The computation of M and J follows the same lines of the one given in detail for the 
low temperature regime. Particularly straightforward is the evaluation of the angular 
momentum, which can be read off from the shift vector The explicit calculation 

gives 


= 



(138) 


where we have used the relation N = L^/2 which holds in our units. Viewing = A 

as a gauge held in 5 dimensions, the angular momentum is equal to the corresponding 
electric charge, as explained in the previous section. The result is then 

7V2 

j = NT-—. (139) 
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This is indeed what we would have expected, since NT is the energy for a gas of N 


particles with Boltzmann density and N‘^/2 is the ground state energy of the N fermions. 


To compute the mass, we used both methods described in sec. 4. Once again, quadratic 
terms in the charge [NT — N‘^/2) appear in the calculation, with different coefficients in 
the two methods. The linear term is however scheme independent and gives the correct 
result 



(140) 


6 Conclusion and open questions 


In this paper we explored the thermodynamic properties of a 1/2 BPS IIB supergravity 
solution called hyperstar. This background was first obtained in |9] by thermal coarse- 
graining of the “bubbling AdS geometry” found in [7]. The hyperstar is in correspondence 
with a distribution of free fermions in thermodynamic equilibrium at temperature T, living 
on a two-dimensional phase space contained in the ten-dimensional geometry. 

We studied both limits of low and high temperature. In the former case, the fermions 
obey the Fermi-Dirac distribution and the supergravity background is obtained from the 
LLM Ansatz by means of a Sommerfeld expansion. We found agreement between the 
energy of the fermions and the ADM mass of the supergravity, modulo a subtlety involving 
terms which we discussed in the main text. We also proposed a way to match the 
entropy of the fermions with the entropy of the hyperstar in the low temperature limit. 
String a' corrections are expected to generate a finite area stretched horizon, lifting the 
naked singularity of the hyperstar to a true black hole singularity. 

In the classical limit of high temperature, we found the explicit form of the metric of 
the supergravity background and we observed how this metric resembles the metric of a 
dilute gas of D3 branes, which corresponds to the S'0(4) invariant sector of the Coulomb 
branch of the CFT. We also computed the associated mass and angular momentum. 

It would be interesting to push this study further. An important point, as already 
remarked, would be to understand better the meaning of the temperature for the super¬ 
gravity solution. On a more fundamental level, it is worthwile to understand the exact 
relation between a thermalized solution like the hyperstar and the Matrix Model descrip¬ 
tion of the half-BPS sector of the dual CFT, extending considerations already made in 


ng. 
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Another issue is whether the appearance of the naked singularity in the hyperstar can 
be understood in terms of a distribution of giant gravitons, as is the case for the superstar 

ng. 

The LLM geometries, upon dimensional reduction to hve dimension, can be seen as 
interesting generalizations of AdS half-BPS extremal black holes [30]. It would then be 
interesting to obtain the explicit dimensional reduction to hve dimension of the hyperstar. 
In this setting one could use the powerful methods of holographic renormalization to carry 
out the computation of the ADM mass. Then one could prove in a rigorous way that the 
quadratic contributions to the mass are effectively spurious and can be eliminated within 
an appropriate renormalization scheme. 

Finally, we would like to mention that the LLM construction has been extended to 
other BPS sectors of type IIB supergravity, see, for instance, [3l] for the 1/4 BPS sector. 
In this case, one modihes the LLM Ansatz in order to accomodate an axion-dilaton held 
which breaks the supersymmetry by half. The ehect of this held is to introduce a dehcit 
angle in the “phase space”. One could try to understand whether this phase space can 
be useful to study the mass and entropy of the corresponding supergravity geometry. 
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